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ABSTRACT 

O The scientific performance of the Planck Low Prequency Instrument (LPI) after one year of in-orbit operation is presented. We describe the main 
^ optical parameters and discuss photometric calibration, white noise sensitivity, and noise properties. A preliminary evaluation of the impact of the 



^ main systematic eff'ects is presented. Por each of the performance parameters, we outline the methods used to obtain them from the flight data and 
^ , provide a comparison with pre-launch ground assessments, which are essentially confirmed in flight. 
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^T) 1 - Introduction bands. The instrument is based on state-of-the-art indium phos- 

, T T 1 . . 1 . 1 1 1 t 1 , phide (InP) cryogenic high electron mobility transistor (HEMT) 

The Planck' mission was designed and developed to produce amplifiers, implemented in a diflferential system using blackbody 

. a full-sky survey of the cosmic microwave background (CMB) ^^^^^ reference signals. The LFI front-end is cooled to ~ 20 K 
^ with unprecedented accuracy, both m temperature and polarisa- ^p^-^^j sensitivity, while the reference loads, connected with 
^ tion. The need to control and remove astrophysical foregrounds fjpj f^ont-end unit, are cooled to ~ 4.5 K. 

I with exquisite precision led to the requirement of intensive spec- n • f^^ u ^ da / ^ 

^— H ^ ^ 1 r ^ • r Following the successiul launch or the Planck satellite 

.. tral coverage over a broad range. P/a^c^ features nine frequency i/iA/r t • ^ u i 

, 1 1 1 . 1 • 11 11 r^rJ. (Kourou, 14 May 2009), the LFI instrument has been work- 

bands, roughly logarithmically spaced, in the range 27-900 GHz. All AU U ^ A A' ^ 

A jj-^- 1 J . .1 • 7 1 -11 -1 mg flawlessly and has been operated according to the mission 

k> As an additional product, the nine Planck sky maps will provide , i n u ^- i \ ^- ^ t onnn 

. . ' . r ^ \ • 1 J . • .1 • plan (Planck Collaboration 201 la). Starting on 4 June 2009, the 

the community with a variety of new astrophysical data in this tct \. a ru ^- j ^ ^- /r^nA7\ 

?-H . , , , y 1 1 / 1 . LFI entered the calibration and performance verification (CPV) 

C3 rich and largely unexplored frequency domain. ^^^.^^^^ instrument testing and radiometer tuning, 

Two complementary cryogemc instruments employing dif- ^ i^jting the different thermal configuration during cool-down 

ferent technologies, the Low Frequency Instrument LFI) and the (Cregorio et al. 201 1). On 13 August 2009, the LFI instrument, 

High frequency Instrument (HFI), share the P/oncfc focal plane ^ ^j^^^ ^^^^ jjPj ^^^^^ ^^j^H-^^ ^^^^^ 

(Bersane h et al. 2010; Lamarre et al. 2010^ The LFI is an ar- ^^^^j observing mode and has been continuously surveying the 
ray of coherent microwave receivers m Ka, Q, and V frequency ^j^^^^ ^^^^ p^^^jl^ pi^^^^^ ^^j^^i^^^ ^^^^^ 

~ ' ~ . , , . . , , ^ egy are given in Planck Collaboration (2011a). 

Planck (http://www.esa.int/P/flncA:) is a proiect of the European xi- -i i ri n-ii 

Space Agency (ESA) with instruments provided by two scientific con- , this paper we provide a description of the in-flight be- 

sortia funded by ESA member states (in particular the lead countries haviour of Planck-LFl and an evaluation of its performance 

France and Italy), with contributions from NASA (USA) and telescope based on the first year of science operations. As planned since 

reflectors provided by a collaboration between ESA and a scientific con- the early design phase, the LFI calibration strategy is based on a 

sortium led and funded by Denmark. combination of ground and in-orbit measurements. Where pos- 
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sible, key instrument parameters have been measured both in 
pre-flight and in-orbit tests, and thus can be compared for self- 
consistency. Ground calibration results have been obtained in 
three main test campaigns, corresponding to three levels of in- 
strument integration, i.e., at single radiometer level (Villa et al. 
2010), at instrument level (Mennella et al. 2010), and during the 
Planck system-level campaign carried out at CSL^. In the lat- 
ter campaign (August 2008), the instruments were operated in 
conditions very similar to flight. 

Planck is based on a novel concept of active cooling, which 
supports for the first time space-borne O.IK bolometers and 
20 K HEMT amplifiers. Much of the complexity of the Planck 
payload comes from its cryo-chain and from the very demand- 
ing stability requirements set by the two instruments (Planck 
Collaboration 201 lb). In particular, monitoring the performance 
and stability of the LFI InP amplifiers and phase switches pro- 
vides ground-breaking technological insight on these devices in 
view of possible future projects. The performance analysis given 
in this paper refers to the first 12 months of data, and represents a 
first assessment of the instrument behaviour and systematics ad- 
equate to support the Planck early release compact source cat- 
alog (ERCSC; Planck Collaboration 2011c) and the early sci- 
ence papers (see Planck papers in this volume). In particular, no 
analysis of polarisation performance is given here. A complete 
characterisation of the LFI instrument behaviour, including a full 
discussion of systematic eflTects and trend analysis, will be pro- 
vided with the January 2013 data release. 

2. Instrument 

2.1. Instrument configuration 

The instrument (Fig. ??) consists of a ~ 20 K focal plane 
unit hosting the corrugated feed horns, orthomode transduc- 
ers (OMTs), and receiver front-end modules (FEMs). A set of 
44 composite waveguides (D'Arcangelo et al. 2009) interfaced 
with the three V-groove radiators (Planck Collaboration 2011b) 
connects the front-end modules to the warm (~ 300 K) back-end 
unit (BEU), which contains further radio frequency amplifica- 
tion, detector diodes, and electronics for data acquisition and 
bias supply. Each LFI radiometer chain assembly (RCA) con- 
sists of two radiometers feeding two diode detectors (Fig. 1), for 
a total of 44 detectors. The 1 1 RCAs are labelled by numbers 
from 18 to 28 as outlined in the right panel of Fig. ??. 

Figure 1 shows a schematic of a complete RCA. The feed 
horn is connected to an OMT, which splits the incoming radia- 
tion into two perpendicular linear polarisation components that 
propagate through two independent pseudo-correlation diflTeren- 
tial radiometers. These radiometers are labelled as M or S de- 
pending on the arm of the OMT they are connected to ("Main" 
or "Side", see lower-left inset of Fig. 1). 

In each radiometer, the sky signal coming from the OMT 
output is continuously compared with a stable 4 K blackbody 
reference load mounted on the external shield of the HFI 4 K box 
(Valenziano et al. 2009). After being summed by a first hybrid 
coupler, the two signals are amplified by ~ 30dB. A phase shift 
alternating at 4096 Hz between 0° and 1 80° is applied in one of 
the two amplification chains. A second hybrid coupler separates 
back the sky and reference load components, which are further 
amplified and detected in the warm BEU. The output voltage 
ranges from -2.5 V to -h2.5 V. 

The output diodes are labelled with binary codes 00, 01 (for 
radiometer M) and 10, 11 (for radiometer S). The four outputs 

^ Centre Spatial de Liege 



of each radiometric chain are referred to as M-00, M-01, S-10, 
S- 11 (Fig. 1). 

After detection, an analog circuit in the data acquisition elec- 
tronics (DAE) removes a programmable oflTset to obtain a nearly 
null DC output voltage, and a programmable gain is applied to 
match the signal level to the ADC input range. 
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Fig. 2. Signal processing after detection with details about the 
analog off'set and gain stages. The analog signal is digitised by 
a 14 bit analog-to-digital converter (ADC) and then processed. 
Processing includes data rebinning, lossy requantisation, loss- 
less compression, and telemetry packetization. 

After the ADC, data are downsampled, requantised, com- 
pressed according to a scheme described in Herreros et al. (2009) 
and Maris et al. (2009), and assembled into telemetry packets. 
On the ground, telemetry packets are converted to volts, using 
the applied off'set and gain factors, and split into sky and refer- 
ence load time-ordered data (TOD). 



2.2. Signal model 

In the ideal case of a perflectly balanced radiometer, the diff'er- 
ential power output for each of the four diodes can be written 
as follows (Seiff'ert et al. 2002; Mennella et al. 2003; BersanelH 
et al. 2010): 

PfuUO = ^ <^tot kfi [Jsky + Tnoise " r (T,,f + Tncse)] , (1) 

where Gtot is the total gain, k is the Boltzmann constant, JS is 
the bandwidth, and a is the diode constant. Tsky and T^^f are the 
sky and reference load antenna temperatures at the inputs of the 
first hybrid, and Tnoise is the receiver noise temperature. In this 
section we always refer to average quantities. We omit, for sim- 
plicity, angled brackets, i.e., T^^y = (T^^y), T^^f = (T^Qf), etc. 
The gain modulation factor, r, is defined by 



r = 



^sky ~^ ^ni 



(2) 



and is used to balance (in software) the temperature off'set be- 
tween the sky and reference load signals and minimise the resid- 
ual l/f noise in the diff'erential datastream. In nominal operating 
conditions, the gain modulation factor is in the range 0.84 < r < 
1.00, depending on channel. This parameter is calculated each 
pointing period^ from the average uncalibrated total power data 
using the relationship: 



r = Vsky/Vref. 



(3) 



^ A pointing period is the amount of time during which the Planck 
telescope scans the same ring in the sky. 
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Fig. 1. Schematic of a complete RCA from feed-horn to analog voltage output, with insets showing the OMT, the details of the 20 K 
pseudo-correlator, and the details of the back-end radio-frequency amplification, low-pass filtering, detection, and DC amplification. 



Details about the gain modulation factor calculation and 
its implementation in the scientific pipeline can be found in 
Mennella et al. (2003) and Zacchei et al. (2011). 

The white noise spectral density at the output of each diode 
is essentially independent of the reference-load absolute temper- 
ature and is given by 



Ar; 



diode 



2 (^sky "I" ^noise) 



(4) 



If the front-end components are not perfectly balanced, then the 
separation of the sky and reference load signals after the second 
hybrid is not perfect and the outputs are mixed. First-order de- 
viations in white noise sensitivity from the ideal behaviour are 
caused mainly by noise temperature and phase- switch amplitude 
mismatches. Following the notation used in Seiff'ert et al. (2002), 
we define ej^ , the imbalance in front end noise temperature, and 
and 6^2, the imbalance in signal attenuation in the two states 
of the phase switch. Equation 4 for the two diodes of a slightly 
imbalanced radiometer then becomes 



^^y.diode^2^^^^diode^2j^^ 



(5) 



which is identical for the two diodes apart from the sign of the 
term {ea^ - ^a^I^-> representing the phase switch amplitude im- 
balance. This indicates that the isolation loss caused by this im- 
balance generates an anticorrelation between the white noise lev- 
els of the single-diode data streams. 

For this reason, the LFI scientific data streams are obtained 
by averaging the voltage outputs from the two diodes in each 
radiometer: 



v; 



rad 



diode 1 



+ W2y; 



diode 2 



(6) 



where w\ and W2 are inverse- variance weights calculated from 
the data as discussed in Zacchei et al. (2011). This way, the 
diode-diode anti-correlation is cancelled, and the radiometer 
white noise becomes 



Ar 



rad 



A 7^ diode 



V2 



(1 + aerj 



1/2 



(7) 



In Eqs. (5) and (7), e «: 1, while a is a term < 1 given by: 



sky 



^ref) 



2 ^^sky "I" ^noise^ (^ref "I" ^noise) 



(8) 



Figure 3 illustrates this anti-correlated noise for a representa- 
tive LFI channel (LFI28M). In black we plot the amplitude spec- 
tral density (ASD) of the weighted sum of the diode pair corre- 
sponding to the radiometer. The sky signal is visible as a series of 
spikes at the satellite rotation frequency and harmonics. The 1// 
and white noise levels are clearly visible, too. In red we show 
the ASD of the weighted diff'erence of the same two diodes, a 
case in which (as expected) the sky signal is removed. In the dif- 
ference signal the white noise is nearly 20% higher because of 
the anti-correlated noise in the output of two diodes. In fact, if 
the diodes were uncorrelated then the sum and diff'erence of their 
ASD's should show identical white noise levels. 

While the eff'ects of this correlated component and the proper 
propagation through the pipeline to maps could be calculated 
and corrected at the map level, we have found it more natural and 
eff'ective to combine the diodes in the time domain, performing 
calibration and further processing on the combined time stream. 

2.3. In-band response and centre frequency 

The in-band receiver response has been thoroughly modelled 
and measured for all the LFI detectors during ground tests. The 
complete set of bandpass curves has been published in Zonca 
et al. (2009). From each curve we have derived the eff'ective cen- 
tre frequency according to 



(9) 



where Ay = Vmax - ^min is the receiver bandwidth and g(y) 
is the bandpass response. Table 1 gives the centre frequencies 
of the 22 LFI radiometers. For each radiometer, g(y) is calcu- 
lated by weight-averaging the bandpass response of the two in- 
dividual diodes with the same weights used to average detector 
timestreams. For simplicity and for historical reasons, we will 



3 



A. Mennella et al.: 




I I J 

0.01 0.01 0.01 0.01 

Frequency (Hz) 

Fig. 3. Amplitude spectral density for weighted sum (black) 
and weighted difference (red) of LFI28M radiometer diodes. A 
white noise component that is anti-correlated between the two 
timestreams accounts for the difference in white noise level seen 
in the plot. Further analysis steps use the weighed sum data. 



Table 1. LFI centre frequencies. 













Radiometer 


M 


Radiometer S 


RCA 


[GHz] 




[GHz] 


V band; "70 GHz" . . . 








LFI18 


71.7 




70.1 


LFI19 


67.5 




69.6 


LFI20 


69.2 




69.5 


LFI21 


70.4 




69.5 


LFI22 


71.5 




72.8 


LFI23 


70.8 




71.3 


Average 




70.3 




Ka band; "44 GHz" . . 








LFI24 


44.4 




44.1 


LFI25 


44.0 




44.1 


LFI26 


43.9 




44.1 


Average 




44.1 




K band; "30 GHz" . . . 








LFI27 


28.3 




28.5 


LFI28 


28.8 




28.2 


Average 




28.5 





continue to refer to the three channels as the 30, 44, and 70 GHz 
channels. 

Colour corrections, C(a), needed to derive the brightness 
temperature of a source with a power-law spectral index a, are 
given in Table 2. The values are averaged for the 1 1 RCAs and 
for the three frequency channels. Details about the definition of 
colour corrections are provided in Zacchei et al. (201 1). 



in flight performance 

3. LFI operations 

3.1. Cooldown and tuning 

LFI operations began during the calibration, performance, and 
verification (CPV) phase of the mission. Functionality and tun- 
ing tests (Gregorio et al. 2011) were carried out, taking advan- 
tage of the varying temperature of the 4K- stage during cooldown 
(see Fig. 4). 



100- 


Commiss. 






LFI Focal plane 
4 K stage 












10 - 




Functionality 










i Pre 
1 tuning . 


Bias tuning 


t 




I 1 1 1 1 1 1 1 


1 1 1 1 1 1 





4Jun2009 9 J un 2009 14Jun2009 19Jun2009 24Jun2009 29Jun2009 

Fig. 4. Cooldown curve of LFI focal plane and 4K stage. 
Functionality and bias tuning tests were carried out during this 
phase. 



For optimal scientific performance, the bias voltages of the 
front-end low noise amplifiers (LNAs) and phase switch currents 
must be carefully tuned. Although radiometer tuning was per- 
formed during ground tests (Cuttaia et al. 2009), the procedure 
was repeated in flight to account for possible changes in the elec- 
trical and thermal environment. 

Phase switch bias currents (/i , /2) of the 30 GHz and 44 GHz 
RCAs were tuned to optimise amplitude balance. This was done 
by exploring a two-dimensional bias surface around the optimal 
points found during ground tests. The results (repeated also at 
the end of the tuning campaign) confirmed the optimal points 
found before launch. The phase switches of the 70 GHz RCAs, 
instead, were set to the maximum biases (1 mA) and were not 
tuned. Ground tests showed that the rise time was sensitive 
to bias currents, and decreased as the current increased from 
0.5 mA to 1 mA (Gregorio et al. 201 1). 

LNA biases were tuned by exploring a large volume in the 
bias space of each amplifier. The common drain voltage, Vd, the 
gate voltage of the first stage, Vgi, and the gate voltage com- 
mon to the remaining stages, Vg2, were sampled according to 
a "hyper matrix" (HYM) tuning strategy, in which several bias 
quadruplets [(Vgi, Vg2)LNAi, (Vgi, yg2)LNA2] were varied for each 
radiometer. This strategy increased considerably the sampled pa- 
rameter space with respect to ground tuning tests (Cuttaia et al. 
2009), and allowed us to fully characterise the radiometer perfor- 
mance in terms of noise temperature, isolation, and drain current 
balance. 

To define the broad bias regions to be deeply sampled by 
the HYM tuning procedure, a pre-tuning phase was run with the 
instrument at 20 K and the reference loads at 20.4 K. The large 
imbalance between the sky and reference load signals provided 
us with enough voltage diflference to estimate noise temperature. 
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Table 2. Colour corrections for the 1 1 LFI RCAs individually and averaged by frequency 



Spectral Index a 



RCA 


-2.00 


-1.00 


0.00 


1.00 


2.00 


3.00 


4.00 


LFI18 


1.054 


1.028 


1.011 


1.003 


1.003 


1.010 


1.026 


LFI19 


1.170 


1.113 


1.066 


1.026 


0.994 


0.969 


0.949 


LFI20 


1.122 


1.079 


1.044 


1.017 


0.997 


0.983 


0.975 


T TTT? 1 


i.Uo/ 




i .UZo 


i .UiU 


1 nnn 
i .uuu 


n QQ^ 


n QQ8 


LFI22 


0.973 


0.971 


0.976 


0.988 


1.007 


1.033 


1.066 


LFI23 


1.015 


1.004 


0.999 


0.998 


1.003 


1.012 


1.026 


70 GHz average .... 


1.070 


1.041 


1.021 


1.007 


1.001 


1.001 


1.007 


LFI24 


1.028 


1.015 


1.007 


1.002 


1.000 


1.003 


1.009 


LFI25 


1.039 


1.024 


1.013 


1.005 


1.000 


0.999 


1.000 


LFI26 


1.050 


1.032 


1.017 


1.007 


1.000 


0.997 


0.997 


44 GHz average .... 


1.039 


1.024 


1.012 


1.004 


1.000 


0.999 


1.002 


LFI27 


1.078 


1.049 


1.026 


1.010 


1.000 


0.996 


0.998 


LFI28 


1.079 


1.049 


1.026 


1.009 


1.000 


0.997 


1.002 


30 GHz average .... 


1.079 


1.049 


1.026 


1.010 


1.000 


0.997 


1.000 



During the HYM tuning phase, smaller bias volumes around 
the optimal pre-tuning points were sampled at four different 
temperatures (see Fig. 4) between 19.1 K and 4K during the 
cooldown of the "^He-JT cooler, allowing us also to character- 
ize the response linearity (Mennella et al. 2009). Drain voltages 
were also tuned for a limited subset of Vgi , Vg2 quadruplets, mak- 
ing the overall bias space six-dimensional. 

Figure 5 shows an example of "condensed noise tempera- 
ture" maps. These are contour plots in the Vgi, Vg2 space for the 
LNAs of a given radiometer. Each point in the plot is the av- 
erage of the best 20% noise temperature values determined by 
the quadruplets sharing that particular Vgi,yg2 pair. The same 
approach was used to map isolation and drain currents. 

Noise temperature and isolation maps were calculated by fit- 
ting data using both a linear and a non-linear response model. 
Minor eff'ects caused by non-linear response were observed, 
as expected, only for 30 and 44 GHz channels, confirming on- 
ground test results. As emphasized in Mennella et al. (2009), 
signal compression in 30 and 44 GHz receivers is relevant only 
if the input range is of the order of > IK. Therefore it can be 
completely neglected in normal operations and in flight calibra- 
tion, where the sky dynamic range is < 10 mK. The optimal 
bias values found during tuning were close to those found dur- 
ing satellite-level tests on the ground, with maximum deviations 
of about 10%. Table 3 sunmiarizes the improvements in noise 
temperature and isolation between the ground and flight tests. 



Table 3. Reduction in Tnoise and increase in isolation (in dB) in 
flight compared to ground tests. 





Min 


Median 


Max 




[%] 


[%] 


[%] 







4.2 


12.5 


J Isolation/Isolation[ JB] . . . . 


1 


6 


115 



The only large change was in radiometer LFI 2 IS, for which 
isolation improved from -7dB measured on ground to -16 dB 
measured in flight. There was a corresponding improvement in 
white noise sensitivity (see Sect. 6.2). 



Table 4. Origin of final bias settings. 



Radiometer 



RCA M S 



LFI 18 Flight^ Flight 

LFI19 Ground^ Ground 

LFI 20 Ground Ground 

LFI21 Ground Flight 

LFI 2 2 Ground Ground 

LFI 2 3 Ground Ground 

LFI24 FEM" FEM 

LFI25 FEM FEM 

LFI26 FEM FEM 

LFI 2 7 Ground Ground 

LFI 2 8 Ground Ground 



^ Determined during in-flight tests. 
^ Determined during satellite-level ground tests. 
Determined during module-level ground tests. 



After tuning, an unexpectedly high level of 1 // noise fluc- 
tuations was observed for the 44 GHz RCAs, using either the 
new flight bias settings or the old ground ones. Dedicated tests 
showed that this instability was correlated with the phase switch 
configuration of 70 GHz LFI23, and disappeared when the ra- 
diometers were biased with the optimal voltages found during 
the optimisation of the individual front-end modules before in- 
strument integration (see Fig. 6). This interaction between RCAs 
belonging to diff'erent frequency channels was unexpected and 
deeply investigated during CPV. The root cause was never fully 
established, but the most likely explanation was a parasitic os- 
cillation triggered by unexpected cross-talk in the warm elec- 
tronics. Details about investigations performed in flight to un- 
derstand and solve this problem are reported in Gregorio et al. 
(2011). The final bias setting (Davis et al. 2009), characterised 
by a slightly higher power consumption (~ 40 mW) but similar 
noise and isolation performance, eliminated the problem. 

Table 4 summarises the bias settings chosen for the front- 
end amplifiers at the end of CPV. They have never been changed 
since the start of nominal operations. 
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Fig. 5. Example of condensed noise temperature maps for the two LNAs of radiometer LFI 2 IS. Contour values represent the ratio 
of the noise temperature in a given bias configuration to the minimum noise temperature found in all tested configurations. The 
yellow crossed point is the chosen bias point. 



10-^ 




Frequency (Hz) 

Fig. 6. Power spectral density of data from LFI2 5M-00 detector 
with two diff'erent bias sets (Blue: optimal biases determined af- 
ter flight tuning tests. Green: optimal biases determined during 
front-end module tuning before integration). 



The last tuning step of the CPV phase configures the signal 
processing unit (SPU) data compressor for each of the 44 LFI 
channels (Maris et al. 2009) so that the data fit into the allo- 
cated telemetry bandwidth without significant loss of sensitivity 
(defined below) from quantisation errors produced by the lossy 
compressor. Quantisation errors appear as an additional white 
noise component that can be characterised as the ratio of the 
quantisation rms Sq over the intrinsic rms cr of the signal before 
compression. The ratio Sqla can be monitored directly in-flight 
with the so-called "calibration channel", a telemetry mode that 
provides 15 min of uncompressed data for each diode every day 
(BersanelH et al. 2010). 

To limit the loss of sensitivity to less than 1%, we require 
Sqicr < 0.1 in the diflTerenced data (Maris et al. 2009). The 



results of the SPU calibration are summarized in Fig. 7 and 
Table 5. The {Sqlo-)^ii < 0.1 requirement in the diflTerenced data 
is satisfied by every channel. 

The stability of the SPU compressor is monitored daily by an 
automatic procedure that checks for variations of the telemetry 
and quantisation error and for near- saturation conditions. Sq has 
changed little during the first year of operations. Typical varia- 
tions are much less than 1%, with maximum variations around 
1.7%. 

3.2. Instrument response during cooldown 

LFI observations began during the cooldown, even before tun- 
ing. Figure 8 shows some of the first diflferential, uncalibrated 
output from LFI28M-00. 

Figure 9 shows a "pseudo-map." The horizontal axis is spin 
axis phase angle, the vertical axis is revolution number, and the 
color scale gives signal amplitude. The map shows the sky, ref- 
erence, and diflTerence data for 100 telescope revolutions on 14 
June 2009. These are the same data that appear in the phase 
binned map in Fig. 8. In the sky signal, one can see the galaxy 
spike, a hint of the CMB dipole, and fluctuations slower than the 
pointing period given by 1 // noise fluctuations in the total power 
voltage output. In the reference load signal, one can see a signif- 
icant gradient due to the rapid cooldown of the 4 K stage during 
this period, as well as the same 1 // fluctuations seen in the sky 
signal. In the difl'erence signal, one can see that correlated fluc- 
tuations in the sky and reference data are cancelled. The mean 
value has been subtracted from each revolution in the diflference 
map. 

3.3. Operations after verification pfiase 

Since the beginning of nominal LFI operations, no instrument 
parameters have been changed, with one exception. Starting on 
11 August 2009, occasional uncommanded jumps in output sig- 
nal in a single channel were seen, sometimes causing saturation 
in the ADC and temporary loss of science data. These jumps 
were traced to single bit-flip changes in the DAE gain- setting 
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LFI18 LFI19 LFI20 LFI21 LFI22 LFI23 LF124 LFI25 LFI26 LFI27 LFI28 

Fig. 7. Average values of Sqla (quantisation error normalized to RMS) measured during the CPV tests for the 44 LFI detectors 
grouped by receiver. For each channel, three bars show Sqlcr for the: (i) sky signal (white); (ii) reference load signal (gray); and (iii) 
differential signal (black). The Sqlcr < 0.1 requirement in the difference data (Maris et al. 2009) is satisfied by every channel. 

Table 5. Normalized quantisation error Sqicr for the N channels at each frequency. Values are plotted in Fig. 7. 



Frequency 










Sqlo- [%] 










Sky 




Reference 




Difference 




Median 


Min 


Max 


Median Min 


Max 


Median Min 


Max 


70 


24 


4.5 


3.1 


5.8 


3.3 2.3 


4.6 


5.1 4.1 


8.0 


44 


12 


10.1 


9.0 


12.4 


12.1 9.2 


15.3 


6.9 4.6 


8.8 


30 


8 


5.7 


5.6 


6.2 


5.1 4.8 


5.3 


5.0 4.6 


5.3 



circuit, presumably a result of single-event upsets (SEUs) from 
cosmic rays. Initially, science data were lost in saturated chan- 
nels until the gain could be reset during the next downlink pe- 
riod. Starting in October 2009, an automatic procedure resets all 
gain values every ~ 40 min. 

Figure 10 shows the cumulative number of gain change 
events as a function of time during the first year of operation. A 
total of 38 were seen through 30 September 2010, correspond- 
ing to an average rate of about one event every 11 days. Out of 
these 38 events, 13 saturated the ADC, leading to lost observa- 
tion time from the corresponding detector, five did not saturate 
the ADC but the compression algorithm was not able to deal 
with the anomalous signal statistics, and 20 could be recovered 
simply by applying the correct gain after the event. 

In most cases the gain increased, but occasionally it de- 
creased. The bit flip occurs in different channels, on different 
bits, and the resulting value after the flip is not always the same. 
The most significant bit (controlling the second amplifier stage) 
never flips, suggesting that the component undergoing a possi- 
ble single event upset (SEU) is the progranmiable gain amplifier 
in the first amplifier stage, an Analog Devices AD526SD/883B. 
The internal technology of the gain amplifiers in the two stages 
is different, therefore a different sensitivity to high energy events 



would not be too surprising. Nevertheless, investigation"^ re- 
vealed no information on SEU effects on AD526 devices. The 
event rate is too low to reveal a correlation between solar ac- 
tivity and cosmic ray fluxes as measured by the onboard space 
radiation environment monitor (SREM; Planck Collaboration 
(2011a)); however, a significant population of high energy cos- 
mic rays, larger than expected before launch, has been observed 
by the HFI instrument (Planck HFI Core Team 201 1). 

3.4. Missing and usable data 

Since the start of the mission, only a small percentage of the 
LFI science data have been lost or considered unusable for sci- 
ence. Table 6 gives the percentage of time lost to missing data, 
anomalies, and maneuvers for the period from 12 August 2009 
to 7 June 2010. 

The main source of missing data is telemetry packets where 
the arithmetic compression performed by the SPU is incorrect, 
causing a decompression error. There were 15 such packets in 
all LFI channels, with negligible scientific impact. For instance, 
for the entire 70 GHz frequency channel, there are just 101 lost 
seconds. 



^ Bibliographic research on Electronic Radiation Response 
Information Center (ERRIC), NASA GFSC Radhome; ESA radi- 
ation effect database 
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Fig. 9. Phase binned uncalibrated pseudo-maps for LFI28M-00. Top: sky signal. Middle: reference load signal. Bottom: difference 
signal. The horizontal axis is spin axis phase angle, the vertical axis is revolution number, and the color scale gives signal amplitude. 
The mean value has been subtracted from each revolution in the difference map. 
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Fig. 8. Phase-binned, uncalibrated, differential data from 
LFI28M-00, acquired on 14 June 2009 near the beginning of 
the performance verification phase. Data for 100 revolutions of 
the telescope have been averaged over 960 "phase bins" of the 
rotation angle. The CMB dipole signal, which was used for pre- 
liminary photometric calibration, and a spike measured while 
crossing the galactic plane, are clearly visible. 
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Fig. 10. Cumulative number of gain change events during the 
first year of operations. 



(2011a)) cause by far the largest fraction of discarded data so 
far; however, we expect these data to be fully recovered after 
additional analysis of the startracker and gyroscope data. If this 
expectation is met, and the performance of Planck remains as it 
has been, well over 99% of all observing time will be usable. 



Anomalies include the DAE gain jumps described in 
Sect. 3.3, and other instabilities (see Sect. 5.1) that make the data 
unsuitable for science. 

Spacecraft maneuvers, from routine repointings of the spin 
axis to stationkeeping maneuvers (see Planck Collaboration 



4. Beams and angular resolution 

The most accurate measurements of the LFI main beams have 
been made with Jupiter, the most powerful unresolved (to 
Planck) celestial source in the LFI frequency range. Since the 
LFI feed horns point to different positions on the sky, they detect 
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Table 6. Percentages of usable and unusable data for the period 
8 August 2009 to 7 June 2010. See text for explanation of the 
categories. 





30 GHz 


44 GHz 


70 GHz 


Category 


[%] 


[%] 


[%] 


Missing 


1.6 X 10-4 


2.7 X 10-4 


3.9 X 10-4 


Anomalies 


0.4 


0.7 


0.4 


Maneuvers 


8.3 


8.3 


8.3 


Used 


91.3 


91.0 


91.3 



Table 7. Dates of Jupiter observations in 2009 



RCA 


Operational Day 


Date 


70 GHz 






18, 23 


168-169 


28 October-29 October 


19, 22 


169-170 


29 October-30 October 


18, 23 


169-171 


29 October-31 October 


44 GHz 






24 


170-171 


30 October-31 October 


25, 26 


163-165 


23 October-25 October 


30 GHz 






27, 28 


170-172 


30 October-1 November 



the signal at different times. Table 7 gives the dates of Jupiter 
observations in the Fall of 2009. 

The first step in extraction of the beams was to remove 
l//-type noise from the data using the Madam destriping map- 
making code. Planets were masked during this process. Details 
of the Madam destriper and of the pipeline implemented to ex- 
tract beams from planet measurements is reported in Zacchei et 
al. (2011). 

To map the beam, each sample contained in the selected 
timelines was projected in the (u, v)-plane perpendicular to the 
nominal line-of-sight (LOS) of the telescope (and at 85° to the 
satellite spin axis). The u and v coordinates are defined in terms 
of the usual spherical coordinates (0, 0): 



u = sin ^ cos 0, 
V = sin ^ sin 0. 



(10) 



To increase the signal-to-noise ratio, data were binned in an 
angular region of 2' for the 70 GHz channels and 4' for the 30 
and 44 GHz channels. We recovered all beams down to -20 dB 
from the peak. An elliptical Gaussian was fit to each beam for 
both M and S radiometers. The FWHM given in Table 8 is the 
square root of the product of the major axis and minor axis 
FWHMs of the individual beams, averaged between M and S ra- 
diometers. The uncertainties in Table 8 are the standard deviation 
of the mean of the Icr statistical uncertainties of the fit. Although 
a small difference between the M and S beams caused by optics 
and receiver non-idealities can be expected, for the purpose of 
point- source extraction relevant for this paper the beams have 
been considered identical, as this effect is well within the statis- 
tical uncertainty. 

The ellipticity and orientation of the beams are fundamental 
parameters. Exhaustive results and details on the determination 
of all LFI beam parameters will be presented in the future. For 
the purpose of this paper. Table 8 gives the typical FWHM and 
ellipticity averaged over each frequency channel. 



Table 8. LFI beam FWHM and mean ellipticity measured in 
flight from the first Jupiter pass. 



FWHM^ Uncertainty'' 
RCA ['] ['] Ellipticity'^ 

70 GHz average 13.01 1.27 

LFI18 13.39 0.170 

LFI19 13.01 0.174 

LFI20 12.75 0.170 

LFI21 12.74 0.156 

LFI22 12.87 0.164 

LFI23 13.27 0.171 

44 GHz average 27.92 1.26 

LFI24 22.98 0.652 

LFI25 30.46 1.075 

LFI26 30.31 1.131 

30 GHz average 32.65 1.38 

LFI27 32.65 1.266 

LFI28 32.66 1.287 

^ The square root of the product of the major axis and minor axis 
FWHMs of the individual RCA beams, averaged between M and S 
radiometers. 

^ The standard deviation of the mean of the Icr statistical uncertainties 
of the fit. Although a small diff'erence between the M and S beams 
caused by optics and receiver non-idealities can be expected, for the 
purpose of point-source extraction relevant for this paper the beams 
are considered identical, as this effect is well within the statistical 
uncertainty. 

^ Ratio of the major and minor axes of the fitted elliptical Gaussian. 



Figure 11 shows three examples of measured beams com- 
pared with calculations performed with the GRASP9^ software. 
The calculated beams have been smeared appropriately to take 
into account satellite rotation and sampling. The eff'ect of sam- 
pling is evident already at -3 dB and cannot be neglected (the 
typical eff'ect on FWHM is 2% at 70 GHz). In the comparison 
of measured and calculated beams, the peak of the simulated 
beam was aligned with the peak of the measured beam (calcu- 
lated from a Gaussian elliptical fit). The electromagnetic model 
of the design telescope (Sandri et al. 2010) was used as a ref- 
erence in the comparison. Ideal parameters have been assumed 
for the shape of the mirrors, the alignment of the telescope and 
focal plane unit, as well as for the pattern of the feed horns. The 
good agreement down to -20 dB demonstrates that, to first or- 
der, the overall LFI optical system is performing as expected. 
Comparison with a more realistic telescope model, including the 
actual alignment, measured mirror shapes, and measured feed- 
horn patterns, will be considered in future work. 

5. Stability and calibration 

5.1. Stability 

Thanks to its differential scheme, the LFI is insensitive to many 
effects in the total-power data caused by 1 // noise, thermal fluc- 
tuations, or electrical instabilities. 

One effect detected during the first survey was the daily 
temperature fluctuation in the back-end unit induced by the 
transponder, which was turned on only for downlinks for the first 
258 days of the mission (Planck Collaboration 2011b, §§2.4.1 



http : //www . ticra . com 



9 



LFI18M 



V 




A. Mennella et al.: LFI in flight performance 



LFI24S 



V 




-o.oas -o.OBts ^0.076 -o.sm) -o.oes -o.oeo 



LFI27M 



V 




-O.OeO -O.OfTS -0.066 -tS.OfiO -O.0S5 



Fig. 11. Example (one for each frequency channel) of the LFI measured beams compared with simulations. The simulated main 
beams have been computed in the co- and cross-polar basis according to the Ludwig's third definition (Ludwig 1973), in spherical 
grids with 301 x 301 points defined with respect to the LOS frame. They are referred to the design telescope configuration. In each 
plot the contours are the levels at -3, -10, -15, and -20 dB from the corresponding power peak. The simulations have been carried 
out in the transmitting mode using GRASP9 software. Physical optics and physical theory of diff'raction has been used on both 
reflectors. 
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Fig. 12. Sky (top), reference (middle) and diflTerence (bottom) 
signals from radiometer LFI27M. The modulation is due to the 
thermal eflTect induced by the satellite transponder being turned 
on and oflf for the daily downlink early in the mission. The quasi- 
sinusoidal -0.1% fluctuation in the sky and reference signals is 
almost completely removed in the diflTerenced data. 



and 5.1). Fig. 12 shows the eflfect of these fluctuations on the ra- 
diometeric output of LFI27M. As expected, the eflfect is highly 
correlated between the sky and reference load signals. In the dif- 
ference, the variation is reduced by a factor ~ (1 - r), where r is 
the gain modulation factor defined in Eq. (3). 



A particular class of signal fluctuations occasionally ob- 
served during operations is represented by electrical instabilities 
that appear as abrupt increases in the measured drain current of 
the front-end amplifiers with a variable relaxation time from few 
seconds to some hundreds of seconds. Typically, these events 
cause a change in both sky and reference load signals and disap- 
pear in the diff'erence (see Fig. 13). 

Because these eflfects are essentially common-mode, their 
residual on the diff'erenced data is negligible, and the data are 
suitable for science production. In a few cases the residual fluc- 
tuation in the diff'erential output was large enough (a few mil- 
likelvin in calibrated antenna temperature units) to be flagged, 
and the data are not used. The total amount of discarded data for 
all LFI channels until Operational Day 389 was about 2000 s per 
detector, or 0.008%. 

A further peculiar eflfect appeared in the 44 GHz detectors, 
where single isolated samples, either on the sky or the reference 
voltage output, were far from the rest. Over a reference period 
of four months, 15 occurrences of single-sample spikes (out of 
24 total anomaly events) were discarded, an insignificant loss of 
data. 



5.2. Calibration 

Photometric calibration, i.e., conversion from voltage to antenna 
temperature, is performed for each radiometer after total power 
data have been cleaned of 1 Hz frequency spikes (see § 7 and 
Zacchei et al. (2011)), and diff'erenced. 

Our calibrator is the well-known dipole signal induced by 
Earth and spacecraft motions with respect to the CMB rest 
frame. The largest calibration uncertainty comes from the pres- 
ence of the Galaxy and of the CMB anisotropics in the measured 
signal. We therefore use an iterative calibration procedure in 
which the dipole is fit and subtracted, producing a sky map that 
is then subtracted from the original data to enhance the dipole 
signal for the next iteration. Typically, convergence is obtained 
after few tens of iterations. 

The thin gray line in Fig. 14 shows the result of this iterative 
process for the radiometer LFI 2 IM, pointing period by pointing 
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Fig. 13. Short spikes in the drain current (left) affect total power signals (right). The jumps are strongly correlated in sky and 
reference signals, so that in the difference data the effect essentially disappears. 



period. In the most stable regions the gain values display relative 
variations of ~ 0.8% rms and ~ 7% peak-to-peak. In the most 
unstable region, where the spacecraft spin was nearly aligned 
with the dipole and the dipole signal was weak, the relative vari- 
ations are ~ 4% r.m.s. and ~ 67% peak-to-peak. These varia- 
tions reflect statistical uncertainties in the determination of the 
gain over a single pointing period, rather than actual changes in 
gain. 

We can put a limit on the true intrinsic gain variations by 
looking at the variation of the total power voltage output (cf. 
§5.1). The small variations in total power constrain intrinsic gain 
variations to be less than 1%. 

The gain solution based on individual pointing periods was 
therefore processed to improve its stability by implementing two 
running averages that have been further smoothed with wavelets: 
one with a 5 -day window, used in the strong dipole regions, and 
the other with a 30-day window, used in the weak dipole re- 
gions. In particular cases, where a true known instrument gain 
change had to be traced^ we used a 5 -day un-smoothed window. 
In Fig. 14 the the thick black line gives the smoothed gain curve; 
variations are now ~ 0.5% rms and ~ 1.6% peak-to-peak. The 
inset in the figures compares the smoothed gain and the relative 
gain variation obtained from Eq. (15), which represent true gain 
fluctuations in the instrument at the level of ±0.2%. In the cur- 
rent gain model implementation these changes are neglected, but 
they will be incuded in future versions of our analysis. 

5.3. Calibration accuracy 

Following COBE and WMAP (Kogut et al. 1996; Jarosik et al. 
201 1), our main calibrator is the dipole modulation in the CMB. 
In our current calibration model we use as calibration signal the 
sum of the solar dipole AJsun and the orbital dipole AT^orb, which 



is the contribution from Planck orbital velocity around the Sun: 



^ e.g., when the satellite started using the transponder always "on" 
which caused a change in the instrument warm unit temperature. 



a: 



(11) 



where ?^axis is the angle between the spacecraft axis and the over- 
all dipole axis (solar -h orbital). 

In Eq. (11), the absolute calibration uncertainty is domi- 
nated by the uncertainty in AJsun, which is known only to about 
1%. The modulation of the orbital dipole by the Earth's motion 
around the Sun is known with an uncertainty almost three orders 
of magnitude smaller; however, at least one complete Planck or- 
bit is needed for its measurement. In the future, calibration based 
on the Earth orbital modulation will be significantly more accu- 
rate. 

The accuracy of our current calibration can be estimated by 
taking into account two components: 1) the statistical uncer- 
tainty in the regions of weak dipole; and 2) the systematic un- 
certainty caused by neglecting gain fluctuations that occur on 
periods shorter than the smoothing window. Neglecting the con- 
tribution from the wavelet filters and assuming, for simplicity, 
that the smoothing is done using plain averages, we can write 
the statistical uncertainty as: 



SG 
~G 



1 n:N{Gi-{G)) 



N-\ 



(12) 



where N is the overall number of pointings over the considered 
period and M is the number of pointings for each smoothing 
region (5 or 30 days). 

Table 9 lists the largest statistical uncertainties in four time 
windows (days 100-140, 280-320, 205-245, 349-389), the first 
two corresponding to minimum and the second two to maximum 
dipole response. 

Our current calibration scheme neglects systematic gain vari- 
ations caused by thermal fluctuations, which introduce spurious 
signal fluctuations on timescales ranging from 1 to 24 hours. 



11 



A. Mennella et al.: LFI in flight performance 




Days after launch 

Fig. 14. Reconstructed gain for radiometer LFI2 IM. Thin gray Hne: gain constants obtained for each pointing period by the iterative 
caHbration procedure. The two regions showing large scatter in the reconstructed gain correspond to regions with minimum ampli- 
tude in the CMB dipole. Thick black line: gain constants obtained after smoothing - this is the actual curve used in the pipeline. The 
inset shows a closer look at the region between Operational Days 200 and 205. Here the relative variation (in %) for the smoothed 
gain model is compared with the relative gain variation calculated from Eq. (15) and shown in dark grey. 



Table 9. Worst-case relative calibration uncertainties. 







[%] 


RCA 


Rad.M 


Rad.S 



70 GHz average 

LFI18 

LFI19 

LFI20 

LFI21 

LFI22 

LFI23 

44 GHz average 

LFI24 

LFI25 

LFI26 

30 GHz average 

LFI27 

LFI28 



0.12 



0.28 


0.24 


0.14 


0.23 


0.25 


0.26 


0.37 


1.04 


0.37 


0.18 


0.29 


0.38 


0.07 




0.35 


0.31 


0.36 


0.30 


0.40 


0.41 



0.05 



0.31 
0.23 



0.35 
0.38 



Being much slower than the satellite spin period, these fluctua- 
tions are well-removed during map-making by the destriping al- 
gorithm (Keihanen et al. 2005, 2010; Zacchei et al. 201 1), so that 
the rms systematic uncertainty per pixel in the final maps due 
to imperfect calibration is at sub-microkelvin levels (see § 7.2). 
However, we expect that a more accurate description of the gain 
variations over short timescales will improve our calibration ac- 
curacy even at the level of time ordered data. If we use the dipole 



signature to estimate an average value Go for the gain over long 
time periods (e.g., several months), we can write the gain versus 
time as 

G(0 = Gox(l+^(0), (13) 

where ^(0 = SG/{G) is the relative gain variation. 

The simplest and most direct measurement of gain changes is 
through the monitoring of the stability of the reference load total 
power voltage output, Vref = K (T^Qf -\- Tnoise)? where K = I /G. 
The relative output voltage variation is: 



SVr.f 
(Vref) 



[(K)j 



(14) 



\(^ref + ^noise)/ \(^ref + ^noise)/ 



1/2 



If the fractional variations in the noise temperature and in the 
reference load signal are negligible compared to gain variations, 
then we have that: 



<Vref> " {K) 



(15) 



Considering that Tnoise is in the range 10-30 K and that the 
reference load temperature is about 4.5 K with fluctuations of 
< 2mK, we have ,^ . < 2 x 10""^. Fluctuations in the 

V ref ■'■^ noise/ 

noise temperature due to thermal changes can be estimated by 
assuming a coupling of 0.5 KK"^ and taking a typical tem- 
perature fluctuation of 1 mK in the front-end unit. We have 
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Fig. 15. Comparison between two gain models for radiometer 
LFI26S. Each model uses the average gain calculated using the 
dipole, but modulates it over time using either Eq. (15) (black 
line) or one of the temperature sensors in the LFI back-end (red 
line). The bottom plot shows the difference between the two gain 
models. 



< 5 X 10 ^, which implies that any variation in the 

(^ref + ^ noise) ^ 

total power signal that is larger than ~ 6 x 10""^ is due to actual 
changes in the gain. 

Figure 15 shows two estimates for ^(t) based on the tem- 
perature sensors in the radiometer back-end and on the relative 
variation of the reference load voltage output (black line). 

In future versions of the LFI gain model, we plan to use the 
iterative solution as a starting point and then trace gain changes 
down to short time scales both by using housekeeping informa- 
tion and by monitoring the relative variation of the total power 
radiometric output (Eq. (15)). 



6. Noise properties 

The noise characteristics of the LFI datastreams are closely re- 
produced by a simple (white -\- 1 //) noise model: 



P(f) = CT^ 



1 -h 



/ 



knee 



(16) 



where P{f) is the power spectrum and a ^ -\. 

In this model, noise properties are characterised by three pa- 
rameters, the white noise limit cr, the knee frequency^ /knee, and 
the exponent of the 1 // component a, also referred to as "slope." 
In this section we show how the LFI noise has been characterised 



^ The frequency at which the 1// and white noise contribute equally 
in power. 
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Fig. 16. White noise level versus time for three radiometers: 
LFI18M, LFI24M, and LFI27M. 



in flight, and we give noise performace estimates based on one 
year of operations. 



6.1. Method 

Noise properties have been calculated following two diff'erent 
and complementary approaches: 1) fitting Eq. (16) to time- 
ordered data for each radiometer; and 2) building normalised 
noise maps by diff'erencing data from the first half of the point- 
ing period with data from the second half of the pointing period 
to remove the sky signal ("jackknife" data sets). 



6.1 .1 . Noise power spectrum estimation in frequency domain 

Estimation of noise power spectra from time ordered data is part 
of the iterative approach in the ROMA map-making suite (Prunet 
et al. 2001; NatoH et al. 2001; de Gasperis et al. 2005). We pro- 
duced joint estimates for both signal (i.e., maps) and noise, and 
then fitted the resulting noise power spectra to Eq. (16) to esti- 
mate the three basic noise parameters cr, /knee, and a. All details 
relative to the pipeline are provided in Zacchei et al. (201 1). 

To verify the stability of the instrument, we produced noise 
estimates from five-day chunks of data separated by 20 days 
each. Figure 16 gives examples for three radiometers. No sig- 
nificant systematic variations were observed in any of the three 
parameters, so that final values and uncertainties could be ob- 
tained by taking the average and the standard deviation of the 
five-day values. 

Results for white noise cr, knee frequency /knee, and slope 
a are reported in Tables 10 and 11. Typical uncertainties are of 
the order of ~ 0.5% for the white noise, between 5 and 10% for 
the slope, and 10 and 20% for the knee frequency. Comparing 
the noise power spectra for all the 22 LFI radiometers with their 
best-fit model, we found good agreement as shown in Fig. 17. 
In the figure, the horizontal red lines represent the white noise 
level (lower line) and the level of equal contribution from white 
and 1 // noise. The frequency corresponding to the intercept of 
the upper red line with the power spectrum is the knee frequecy 

/mee- 
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Fig. 17. Power spectral densities (PSDs) for the LFI radiometers. The lower red line marks the white noise level. The upper red line 
marks the level of equal contribution from white and 1 // noise, so the intercept with the power spectrum marks the knee frequecy, 
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Table 11. Knee frequency and slope for the 22 LFI radiometers. 

Knee Frequency Slope 

Radiometer M Radiometer S Radiometer M Radiometer S 
RCA [mHz] [mHz] 



70 GHz 



LFI18 


16.3 + 3.9 


LF119 


15.1 + 2.2 


LF120 


18.7+ 5.3 


LFI21 


37.2+ 6.1 


LFI22 


12.7 +19.7 


LF123 


34.6 + 4.8 



44 GHz 



LFI24 


46.2 + 


3.9 


LFI25 


24.9 + 


3.7 


LFI26 


67.6 + 


2.2 



30 GHz 

LFI27 187.4 +29.5 

LFI28 122.2+10.2 



Table 10. White noise sensitivities for the 22 LFI radiometers. 





White Noise Sensitivity 




Radiometer M 


Radiometer S 


RCA 


[yUKcMBsl/2] 


[yL/KcMBS^/^] 


70 GHz 






LFI18 


512.0+1.8 


465.7+1.7 


LFI19 


581.4+1.6 


555.6+2.4 


LFI20 


590.8+1.7 


623.2+1.5 


LFI21 


455.2 +2.8 


564.1+6.2 


LFI22 


492.0+2.3 


534.4 +2.8 


LFI23 


507.7+1.7 


542.4 ±2.6 


44 GHz 






LFI24 


462.2+1.9 


399.2+1.3 


LFI25 


413.6+3.6 


392.6+1.6 


LFI26 


478.6+3.1 


418.6+4.2 


30 GHz 






LFI27 


277.7 +2.1 


302.9+1.6 


LFI28 


312.3+1.7 


285.3+1.4 



6.1 .2. Jackknife noise maps 

During each pointing period Planck repeatedly scans essentially 
the same circle in the sky. This provides a powerful method to 
remove the sky signal and produce maps containing only the in- 
strumental noise (both uncorrected and correlated on timescales 
shorter than about 20 minutes). The details of this method are 
given in Zacchei et al. (201 1). These "jackknife" noise maps can 
be normalised to the white noise estimate at each pixel obtained 
from the white noise covariance matrix, so that a perfectly white 
noise map would be Gaussian and isotropic with unit variance. 
Figure 1 8 shows an example of such normalised noise maps for 
the 30 GHz LFI frequency channel. The figure shows a generally 
structureless map, as expected, apart from a few regions in the 
Galactic plane where the temperature maps have large gradients 
over the pixel scale, causing the sky signal to leak into the noise 
map. In subsequent analysis steps we applied a ~ 20% mask 
to remove these regions (see Fig. 19). Table 12 gives the stan- 
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-1.16 + 0.01 
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104.4 +26.1 
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-0.91 +0.13 




Fig. 19. 30 GHz mask used to extract noise parameters from the 
corresponding noise map (masked points are shown in blue). 

Table 12. Standard deviation of normalised noise maps. 



30 GHz 


44 GHz 


70 GHz 


1.039 


1.016 


1.002 



dard deviation of normalised noise maps obtained for the three 
LFI frequencies. Deviations from unity trace the contribution of 
residual 1// noise in the final maps, which ranges from 0.2% at 
70 GHz to 4% at 30 GHz. A further Gaussianity test performed 
by calculating skewness and kurtosis of the normalised map af- 
ter masking the galactic plane yielded null values within two 
standard deviations of 500 Gaussian simulations. 

6.2. Sensitivity 

Fig. 20 compares the calibrated white noise sensitivity for the 22 
LFI radiometers calculated from flight data and during ground 
tests performed both at instrument and satellite level. Error bars 
reported for values measured in flight and during satellite tests 
are statistical uncertainties derived from calculations performed 
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-4.0 4.0 Sigma 

Fig. 18. 30 GHz normalised jackknife noise map. We also highlight a region in the Galactic plane affected by leakage from strong 
foreground emission. The inset at high galactic latitude highlights the white nature of the noise. 



on several data chunks^. Error bars reported for values measured 
at instrument level, instead, represent the uncertainty between 
two different methods used to extrapolate the white noise sen- 
sitivity measured at test conditions (~ 20 K input and ~ 26 K 
front-end temperature) to flight conditions (~ 2K input and 
~ 20 K front-end temperature). Further details about this extrap- 
olation are reported in Mennella et al. (2010). 

Figure 20 shows general agreement in the sensitivity cal- 
culated in various test campaigns, with two outstanding excep- 
tions, radiometers LFI24M and LFI21S, which displayed signif- 
icantly improved noise levels in-flight. This resulted from an in- 
correct bias setting of these radiometers during ground tests (see 
Mennella et al. 2010) and was resolved during in-flight tuning. 

In Fig. 20 values for LFI18M and LFI24M measured at 
instrument-level are not reported because these two radiometers 
failed. LFI18M was replaced with a spare unit and LFI24M was 
repaired before instrument delivery (see Mennella et al. 2010). 

Table 13 summarizes the sensitivity numbers calculated dur- 
ing the first year of operations using methods and procedures 
outlined in § 6.1 and described in detail in Zacchei et al. (2011) 
compared with scientific requirements. The measured sensitiv- 
ity in very good agreement with pre-launch expectations. While 
the white noise moderately exceeds the design specification, this 
performance is fully in line with the LFI science objectives. 

The consistency of the noise estimates outlined in Table 10 
was tested by comparing the jackknife noise maps, the white 
noise covariance matrices, and 101 noise Monte Carlo realisa- 
tions (see Zacchei et al. 2011, for details). In fact, noise covari- 
ance matrices and Monte Carlo simulations are dependent on 
estimated noise parameters, while the jackknife noise maps de- 
scribe directly the noise in the final maps. 



^ About 20 one-hour chunks for on-ground satellite data, one year of 
operations for in-flight data. 



Table 13. White noise sensitivities for the LFI frequency chan- 
nels compared with requirements. 



Channel 


Value 


Requir. 




70 GHz 


152.6^KcMB s^/' 


119^KcMB 


,1/2 


44 GHz 


173.1^Kcmbs1/2 


119a/Kcmb 


,1/2 


30 GHz 


146.8^Kcmbs1/2 


119yL/KcMB 


,1/2 



We produced pseudo-Q spectra from the jackknife noise 
maps and noise Monte Carlo maps by anafast (see Fig. 21), 
and compared high multipole tails (1 150 < ^ < 1800) to the pre- 
dictions from the white noise covariance matrices. In the analy- 
sis we masked out a 20% region of the galactic plane and un- 
solved pixels (pixels that have HEALPix^ bad pixel value in 
noise maps), leaving us with sky fraction /sky 0.8. 

Figure 22 shows a comparison of the noise estimates ob- 
tained with the various methods. In the figure the green trian- 
gles refer to noise obtained from binned pure white noise Monte 
Carlo maps, i.e., maps containing no residual 1// noise. At all 
frequencies the high-^ mean of the jackknife noise map pseudo- 
Q spectra weighted by the inverse sky coverage, (/sky)~^ (shown 
in red), are in the 68% range of the noise Monte Carlos (rep- 
resented with the back error bars). The point from the 30 GHz 
jackknife noise map is almost l-cr higher than the median value 
for the noise MC due to residual gradient leakage (due to 
point sources etc.), which is strongest for the 30 GHz channel. 
The noise power from the white noise covariance matrices and 
binned maps is always lower than the noise power from the full 
noise Monte Carlo and jackknife noise maps, as expected, due 

^ http://healpix.jpl.nasa.gov 
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Fig. 20. Comparison of white noise sensitivity values calculated in flight and on ground during the satellite-level and instrument- 
level test campaigns. Values for LFI18M and LFI24M measured at instrument-level are not reported because these two radiometers 
failed before instrument delivery and were replaced and repaired, respectively. 
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Fig. 21. The noise angular power spectra. Black: noise Monte Carlo (median, 16% quantile, and 84% quantile). Red: jackknife noise 
map pseudo-Q. Blue: values from white noise covariance matrices. 



to residual 1 // noise present even at high multipoles. This eflPect 
is especially pronounced in the 30 GHz channels, which have 
larger knee frequencies (see § 6.3). 



6.3. 1/f Noise properties 

Figure 23 shows a comparison of the knee frequency for the 
22 LFI radiometers estimated from flight data and ground satel- 
lite calibration. Knee frequency and slope after one year of op- 
erations for all LFI radiometers were given earlier in Table 1 1 . 
With some exceptions the values are below or compatible with 
requirements and are comparable within the error bars among 
diflTerent measurements. 

There are a few cases (LFI23S, LFI24S, LFI27M, LFI28M, 
LFI27S) in which the knee frequency measured in flight is 
higher than that measured on ground. The cause is still under in- 
vestigation, but the impact on the scientific quality of the data is 
small. Destriping eff'ectively removes correlated structures, lim- 



iting the impact of the high knee frequencies to a small (< 4%) 
increase of the noise variance (see Table 12 and Fig. 18). 

As shown in the low-^ part of the angular power spectra in 
Fig. 21, the noise parameters estimated in flight are a good rep- 
resentation of the noise in the actual maps. This means that for 
cosmological purposes we cannot rely on simple white noise co- 
variance, and a full timeline-to-map Monte Carlo is needed. On 
the other hand, this also means that there is no need for a more 
complex functional noise model with respect to the one reported 
in Eq. (16). 



7, First assessment of systematic effects 

In this section we discuss the most relevant systematic efl'ects 
in the first year data. The LFI design was driven by the need to 
suppress systematic eflfects well below instrument white noise. 
The receiver diff'erential scheme, with reference loads cooled 
to 4 K, greatly minimises the eff'ect of 1 // noise and common- 
mode fluctuations, such as thermal perturbations in the 20 K LFI 
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Fig. 23. Comparison of 1 // knee frequencies measured in flight and on ground during the satellite-level test campaign. 



focal plane. The use of a gain modulation factor (see Eq. (3)) 
largely compensates for spurious contributions from input oflT- 
sets. Furthermore, diode averaging (Eq. (6)) allows us to cancel 
second-order correlations such as those originating from phase 
switch imbalances. 

We have developed an error budget for systematic eflTects 
(Bersanelli et al. 2010) as a reference for both instrument de- 
sign and data analysis. Our goal is to ensure that each systematic 
eff'ect is rejected to the specified level, either by design or by 
robust removal in software. At this stage, the following eflTects 
proved to be relevant: 

- 1// noise (discussed in § 6.3), 

- 1 Hz frequency spikes, 

- thermal fluctuations in the back-end modules driven by tem- 
perature oscillations from the transponder during the first 
survey, 

- thermal fluctuations in the 20 K focal plane, 

- thermal fluctuations of the 4 K reference loads. 

For each of these eff'ects we used flight data and information 
from ground tests to build timelines, maps, and angular power 
spectra that represent our best knowledge of their impact on the 
scientific analysis. Figure 24 shows maps and histograms of the 
combination of all the above mentioned systematic eff'ects for 
the three LFI channels. The figure shows that the 30 and 44 GHz 
channels display a residual level of spurious signals about twice 
larger than the 70 GHz channel. As shown in the following sec- 
tions, this is caused by temperature fluctuations of 4K refer- 
ence loads that are significantly larger at 30 and 44 GHz than 
at 70 GHz. 

A convenient way to assess the eflfect on angular power spec- 
tra is to calculate the ratio = Cf^^/C^^^^^, where Cf^^ is the 
angular power spectrum of the systematic eflfect map and C^^^^^ is 
the angular power spectrum of the instrument noise contribution 
including the residual 1 // component remaining after destrip- 
ing. Figure 25 shows for the three LFI frequency channels. 
In each panel we have plotted the spectrum obtained from the 
global maps in Fig. 24 and the spectra derived from the single- 
component maps (see next sections). 



Table 14. Summary of residual eff'ects on maps in //Kcmb from 
main systematic eff'ects. 



30 GHz 44 GHz 70 GHz 





P-P 


r.m.s. 


P-P 


r.m.s. 


P-P 


r.m.s. 


1-Hz spikes 


4.00 


0.45 


1.51 


0.15 


2.56 


0.30 


BEM T fluct. 


1.27 


0.11 


0.63 


0.05 


2.70 


0.24 


FEM T fluct. 


1.05 


0.23 


1.15 


0.22 


1.12 


0.21 


4KT fluct 


9.76 


0.98 


9.73 


0.98 


1.30 


0.16 


Total^ 


10.92 


1.10 


9.73 


0.98 


4.28 


0.45 



The total has been estimated from the maps combining all systematic 
eff^ects (Fig. 24). 



Figure 25 shows that is in the range 10" ^-10"^ for £ in the 
range 2-20, and p^ < 10"^ for ^ > 10. The figure also shows that 
at 30 and 44 GHz the largest contribution to systematic eflfects 
is determined by temperature fluctuations of the 4 K reference 
loads, while the residual systmatic uncertainty in the 70 GHz 
channel is mainly caused by back-end temperature fluctuations 
and, at small angular scales, by frequency spikes. Table 14 gives 
an overview of our current assessment of residual peak-to-peak 
and rms systematic eflTects per pixel on LFI temperature maps. 
Maps were made with A^side = 5 12 at 30 GHz and N^idc = 1024 at 
44 and 70 GHz. Corresponding pixel sizes are ~ 6.'8 and ~ 3. '4. 

Further advances in the data analysis pipeline will be aimed 
at removing spikes from the 70 GHz data, improving relative cal- 
ibration to account for thermally driven fluctuations, and further 
suppressing spurious fluctuations caused by 4 K temperature in- 
stabilities at 30 and 44 GHz. 

7.1. Frequency spikes 

Spikes are seen in the radiometer outputs in the frequency do- 
main at multiples of 1 Hz. These "frequency spikes" were first 
detected during ground tests, and are caused by pickup from 
the clock of the housekeeping electronics (Meinhold et al. 2009; 
Mennella et al. 2010). The pickup occurs between the detector 
diodes and the DAE gain stage. Frequency spikes are present at 
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Fig. 24. Maps of the combined systematic effects from frequency spikes and thermal fluctuations. 



some level in the output from all detectors, but aflTect the 44 GHz 
data most strongly because of the low voltage output and high 
DAE gain values in that channel. In this section we provide es- 
timates of the residual systematic eflTect after frequency spike 
removal at 44 GHz, and the eflTect caused by the spikes without 
any removal at 30 and 70 GHz. 

In the time domain, the frequency domain spikes comprise 
a one second square wave with a rising edge near 0.5 s and a 
falling edge near 0.75 s in on-board time. During the first year of 



operations we did not observe any deviation in either the phase 
or the shape of these signals. 

Templates of these spurious square waves obtained from the 
output voltages have been used to remove this eflTect from the 
data before diflferentiation (procedures and algorithms are de- 
scribed in Zacchei et al. (201 1)). Although we could in principle 
apply the removal process to all data, we decided to remove the 
spikes only from the 44 GHz data, which are the most aflfected 
by this eflTect. 
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Fig. 22. Comparison of noise estimates from different methods. 
Black: noise from high-^ means of noise Monte Carlo statistics. 
Red: noise from the high-^ mean of the jackknife noise map 
pseudo-Q spectra. Green: noise from the binned white noise 
Monte Carlo maps. Blue: noise from white noise covariance ma- 
trices. Error bars are relative to values from Monte Carlo simu- 
lations. We have weighted all the values with the inverse of the 
analysed sky fraction, f~^, to represent the noise levels of tem- 
perature Q spectra that would result from using full sky maps. 



To estimate the residual in the spike removal process caused 
by slow variations in the square wave amplitude, we used a sim- 
ple minimisation procedure to estimate the amplitude of the 
template for each hour of data, then smoothed the amplitude 
by simple binning with a 10-day window size. We took this 
smoothed amplitude as our estimate of the true time variation 
of the spike signal, and calculated the residual error that would 
result if we approximated this time-varying signal with a con- 
stant template in the spike removal process. Although we could 
use, in principle, the time-varying template to remove the effect, 
we decided to use the simplest and most robust approach of a 
constant-amplitude template. The possibility to implement of a 
model that accounts for slow spike amplitude drifts will be con- 
sidered for future data releases. 

Templates were also used to generate maps for each fre- 
quency channel to determine the impact of the spike signal on the 
science results. These maps and the corresponding histograms 
are shown in Fig. 26. For the 44 GHz channel we also show 
the residual effect after removing the constant template, show- 
ing that after removal the residual effect is reduced by a factor of 
about 20. 




Fig. 25. Ratio = Cf^^/C^^^^^ for the three LFI frequency chan- 
nels (top: 30 GHz; middle: 44 GHz; bottom: 70 GHz). In each 
panel we have plotted both the global spectrum obtained from 
the maps in Fig. 24 and the individual spectra calculated from 
the single-component maps. 
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Fig. 26. Simulated spike maps for each channel, and the residual after removal for the 44 GHz channel. Maps have been obtained 
using constant- amplitude templates to generate spike TODs. The datastream used to generate the residual map has been obtained 
by subtracting a spike TOD with time- varying amplitude from a constant-amplitude spike TOD. 
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7.2. Thermal fluctuations 

The LFI is sensitive to temperature fluctuations of the warm 
back-end unit, the 20 K focal plane, and the 4 K reference loads. 
In the first two, temperature variations impact the sky and ref- 
erence load signals at a similar level, so that in the diflTerential 
radiometric output the residual spurious variations are reduced 
by more than one order of magnitude. Fluctuations at the level 
of the 4 K reference loads, instead, transfer directly to the ra- 
diometric output, and may represent a more critical source of 
systematic errors. 

The eflTects of low frequency thermal fluctuations are 
strongly suppressed by the spacecraft spin itself and the de- 
striping map-making codes (Zacchei et al. 2011). High fre- 
quency thermal fluctuations are strongly damped by the thermo- 
mechanical structure of the spacecraft and instruments (Tomasi 
et al. 2009). The dominant efl'ects, then, come from fluctuations 
at frequencies in a range around the spin frequency. 

In this section we present an overview of the impact on the 
LFI science caused by the residual level of systematic efl'ects in 
the final maps using a dedicated pipeline that follows the proce- 
dure outlined below. 

- Start from a time ordered datastream of a temperature sensor 
representative of the temperature behaviour of the consid- 
ered thermal stage; 

- low-pass filter to remove high-frequency sensor noise; 

- apply thermal transfer functions where appropriate to ob- 
tain the physical temperature behaviour at the level of the 
receiver components sensitive to the fluctuation; 

- apply the radiometric transfer function to convert the physi- 
cal temperature fluctuation into antenna temperature fluctu- 
ation; 

- resample at scientific sampling frequency, calibrate and build 
diflferenced time-ordered data; 

- build maps with Madam. 

These maps have been generated combining actual flight 
housekeeping data with thermal and radiometric transfer func- 
tions obtained from flight and ground- test data (Terenzi et al. 
2009). They represent, therefore, our current best estimate of the 
impact of the individual efl'ects on the science in flight condi- 
tions. 

As explained in § 5.3 further developments of this work will 
aim at including thermal fluctuations from the back-end unit in 
the calibration model which will reduce the need to remove the 
spurious signal from the time ordered data during map-making. 

7.2.1 . Back-end temperature fluctuations 

As mentioned in § 5.1, for the first 258 days of the mission the 
satellite transponder was switched on only for the daily telecom- 
munication period. This induced quasi- sinusoidal fluctuations in 
the Planck service module temperature that propagated to the 
LFI back-end unit, and were recorded by the instrument house- 
keeping sensors at the level of ±0.2 K/day. 

This temperature variation drives variation in the total power 
voltage output of both sky and reference detectors, which largely 
disappears in the difl'erence, shown in Fig. 12. The efl'ect on the 
undiflferenced total power data, however, was useful in calcu- 
lating correlation coeflftcients between temperature and voltage 
outputs. These coefl&cients have been used to estimate the resid- 
ual efl'ect on time-ordered data and maps. Figure 27 shows this 
correlation for the detector LFI27M-00. 



1.45 
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Fig. 27. Correlation between back-end temperature and total 
power voltage output relative to the LFI detector LFI27M-00. 



Figure 28 shows the expected systematic efl'ects of temper- 
ature fluctuations of the back-end unit during the first year of 
operations. The peak-to-peak eff'ect is ~ 1 yuKcMB, and the rms 
is well below 1 //Kcmb • 

7.2.2. Front-end temperature fluctuations 

Another source of temperature fluctuations is the Planck sorp- 
tion cooler (Morgante et al. 2009; Planck Collaboration 2011b), 
which displays temperature variations of the order of ~ 0.5 K 
peak-to-peak at the cold end attached to the LFI focal plane unit. 
These fluctuations are damped by a thermal stabilisation assem- 
bly (TSA; Planck Collaboration 2011b) down to about 100 mK 
peak-to-peak at low frequency. 

During the second half of the first year of operations the ex- 
pected cooler degradation led to an increase of the cold end tem- 
perature and the low-frequency fluctuations. At the radiometers 
the temperature increased by about 0.2 K, with fluctuations run- 
ning about 30 mK peak-to-peak. We estimated the efl'ect of these 
fluctuations on the science data using ground-measured trans- 
fer functions that include the eff'ect of temperature on both gain 
and noise. Figure 29 shows the calculated efl'ects for the three 
LFI frequency channels. As in the case of back-end temperature 
fluctuations, the overall eff'ect is ~ 1 //Kcmb peak-to-peak and 
0.2yuKcMB rms. 

7.2.3. 4K reference load temperature fluctuations 

The temperature stability of the 4 K reference loads, attached to 
the HFI 4 K mechanical box, is a key factor in the LFI systematic 
eflfects budget. Temperature fluctuations of the loads are driven 
primarily by fluctuations at the 20 K cold-end interface, which 
propagate to the LFI reference loads by conduction through the 
HFI mounting struts and radiation from the LFI. See Planck 
Collaboration (2011b) for details. 

The 70 GHz reference loads are mounted close to the HFI 
4 K plate and are stabilised by proportional-integral-derivative 
(PID) active controllers that remove almost completely low fre- 
quency fluctuations. The loads at 30 and 44 GHz are farther from 
the HFI 4K plate and closer to the HFI mounting struts; their 
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Fig. 28. Simulation of the systematic effects of temperature fluctuations of the LFI warm back-end unit. Maps are A^side = 512 at 
30 GHz, and A^side = 1024 at 44 and 70 GHz. 



temperature is less stable. Details of temperature sensor mea- 
surements are provided in Planck Collaboration (2011b). 

This diflTerence is reflected in the maps shown in Fig. 30, 
where the residual systematic eflTect in the 30 and 44 GHz 
channels is ~ 10 times larger than in the 70 GHz channel. 
Nevertheless, given the reduction provided by destriping, the 
residual systematic is several times below the instrument noise, 
as shown by the angular power spectra in Fig. 25. Further work 



will be aimed at correlating radiometric and housekeeping data 
to further reduce the residual eflTect caused by these fluctuations. 

8. Summary of main performance parameters 

Table 15 gives a top-level summary of instrument performance 
parameters measured in flight. Optical properties have been suc- 
cessfully reconstructed using Jupiter transits and the main pa- 
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Fig. 29. Simulation of the systematic effects of temperature fluctuations of the LFI focal plane unit. Maps are A^side = 512 at 30 GHz, 
and A^side = 1024 at 44 and 70 GHz. 



rameters are in agreement with pre-launch estimates. Also the 
white noise sensitivity agrees with ground measurements; for 
some channels sensitivity improved after in-flight bias tuning. 
Parameters describing the 1 // noise component are in line with 
ground measurements and the 50mHz requirement except at 
30 GHz. That channel that has a knee frequency over 100 mHz, 
which is, however, well-handled by destriping. Absolute photo- 
metric calibration based on the CMB dipole yields an overall 



statistical uncertainty of < 1%. Variations due to slow instru- 
mental variations are traced by the calibration pipeline, yield- 
ing an overall uncertainty ranging between 0.05% and 1%. The 
residual systematic uncertainty is of the order of 1 yuK rms per 
pixel. Average colour corrections are provided in Table 2. 
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Fig. 30. Simulation of the systematic effects of temperature fluctuations of the 4 K reference loads. Maps are N^ide = 512 at 30 GHz, 
and A^side = 1024 at 44 and 70 GHz. Note that the scale of 30 and 44 GHz maps is 10 times larger than then scale of the map at 
70 GHz 



9. Conclusions 

We have discussed the scientific performance of the Planck Low 
Frequency Instrument after one year of operations. 

Since the start of Planck nominal operations, the LFI has 
shown excellent stability in all measured parameters. The instru- 
ment uninterruptedly observed the microwave sky with negligi- 
ble data loss and less than 1 % discarded data because of anoma- 



lies. Typical variations in the instrumental output were less than 
1% on time scales of several days and were mainly driven by 
slow thermal fluctuations. 

The main beams have been reconstructed down to -20 dB 
using Jupiter as a source, with results closely matching those ex- 
pected from simulations. In-flight measurements of white noise 
sensitivity are in very good agreement with ground results, with 
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Table 15. Summary of main LFI performance parameters. 



Parameter 


JULrrlZ 


44 LrrlZ 


/U LrrlZ 


Centre frequency [GHz] 


28.5 


44.1 


70.3 


FWHM ['] 


32.65 


27.92 


13.01 


Ellipticity 


1.38 


1.28 


1.26 


Whifp, Tini<;p, <5P,Ti<jitivitv r//PC/^**T3 <5^/^l 


146 8 


173.1 


152.6 


/knee [mHz] 


113.7 


56.2 


29.5 


1// slope 


-0.87 


-0.89 


-1.03 


Absolute calibration uncertainty[%] 


1 


1 


1 


Relative calibration uncertainty [%] 


0.05 


0.07 


0.12 


Systematic eff'ects uncertainty [//Kcmb] 


1.10 


0.98 


0.45 



significant improvements in some channels thanks to in-flight 
bias tuning. The impact of low frequency noise is very small, es- 
pecially at 70 GHz where the measured knee frequency is about 
30 mHz, almost twice better than the requirement. At 30 GHz 
the 1// noise is higher than measured on the ground, with 
a knee frequency of about 100 mHz. Residual fluctuations in 
the timestreams, however, are eff'ectively removed during map- 
making, as expected from several pre-launch simulation studies 
(Kurki-Suonio et al. 2009; Poutanen et al. 2004; Keihanen et al. 
2004). We have shown that the noise increase due to the l/f 
component is a few percent at 30 and 44 GHz, and only 0.2% at 
70 GHz. 

Photometric calibration is based on the CMB dipole via an 
iterative procedure explained in Zacchei et al. (2011). Excellent 
absolute calibration accuracy of < 1% will likely be improved 
still further in future analyses by using the orbital dipole as abso- 
lute calibrator. Our current relative calibration traces gain varia- 
tions on timescales larger than 5-10 days, yielding an overall sta- 
tistical accuracy in the range 0.05-0.1%. Thermally-driven gain 
fluctuations on smaller timescales are currently not implemented 
in our gain model and contribute as a systematic uncertainty in 
the final maps. A new version of the gain model, now being de- 
veloped, will take into account the effect of such fluctuations to 
further reduce this residual uncertainty. 

We have presented a preliminary analysis of all the system- 
atic effects that are relevant at this stage of the analysis. Their 
impact on LFI science has been evaluated by projecting each ef- 
fect on full-sky maps and angular power spectra through a ded- 
icated pipeline that exploits radiometric and housekeeping data 
in conjuction with external parameters measured in ground or 
flight tests, such as thermal susceptibilities or spike character- 
istics. The combined residual effect from frequency spikes and 
thermal instabilities (in the 4 K, 20 K, and 300 K stages) is sev- 
eral order of magnitudes below the instrumental noise at all an- 
gular scales smaller than 10°, and less than 10% at larger scales. 

The overall performance of LFI as measured in-flight 
demonstrates an instrument with unprecedented combination of 
sensitivity, angular resolution and suppression of systematic er- 
rors for full- sky imaging at these frequencies. In particular, the 
control of systematic effects and non-white noise components 
represent key challenges for accurate extraction of the cosmo- 
logical information encoded in the temperature and polarisation 
maps. Our preliminary assessment shows that, even without ded- 
icated de-correlation of thermal effects, the LFI is already largely 
immune to instrumental effects, with prospects of further sup- 
pression after implementing a more representative gain model 
and temperature fluctuation removal algorithms. 



One more year of continuous observations is currently 
planned for Planck, with a further LFI extension to the maximum 
lifetime allowed by the sorption cooler. Everything to date sug- 
gests that the instrument will maintain its performance through- 
out the remaining period and provide rich and high-quality sci- 
entific data that will be explored for many years to come. 
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